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Abstract. We study the conditions for convective instability in rotating, non-magnetic proto-neutron stars. The criteria that 
determine stability of nascent neutron stars are analogous to the Solberg~H0iland conditions but including the presence of 
lepton gradients. Our results show that, for standard angular velocity profiles, convectively unstable modes with wave-vectors 
parallel to the rotation axis are suppressed by a stable angular momentum profile, while unstable modes with wave-vectors 
perpendicular to the axis remain unaltered. Since the wave-vector is perpendicular to the velocity perturbation, the directional 
selection of the unstable modes may result in fluid motions along the direction of the rotation axis. This occurs in rigidly 
rotating stars as well as in the inner core of differentially rotating stars. Our results provide a natural source of asymmetry for 
proto-neutron stars with the only requirement that angular velocities be of the order of the convective characteristic frequency. 
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1. Introduction 



Hydrodynamical instabilities in newly born neutron stars play 
an important role in enhancing neutrino luminosities and in- 
creasing the energy deposition efficiency. It was first argued by 
Epstein (1979) that the negative lepton gradient that arises in 
the outer layers of the proto-neutron star (PNS) can be con- 
vectively unstable. However, convection in PNSs can be driven 
not only by the lepton gradient but by a negative entropy gradi- 
ent as well (see, e.g., Arnett 1987, Burrows & Lattimer 1988). 
The development of negative entropy and lepton gradients is 
common in many simulations of supernova models (Bruenn & 
Mezzacappa 1994, Bruenn et al. 1995) and evolutionary mod- 
els of PNSs (Burrows & Lattimer 1986, Keil & Janka 1995, 
Sumiyoshi et al. 1995, Pons et al. 1999). The nature of in- 
stabilities in PNSs has been considered by a number of au- 
thors (Grossman et al. 1993; Bruenn & Dineva 1996; Miralles 
et al. |200Q| MPU in the following; Miralles et al. 120021 In 
non-rotating and non-magnetic PNSs, the criteria indicate the 
presence of two essentially different instabilities (MPU) with 
a convectively unstable region surrounded by a neutron-finger 
unstable region, the latter involving typically a larger portion of 
the stellar material. The unstable zones grow in size and move 
inward as the PNS deleptonizes (Keil et al. 1996), until neutrino 
diffusion progressively reduces the temperature and lepton gra- 
dients, and convective instabilities disappear, at most in 30-40 
s. The growth time of convective and neutron finger instabilities 
are substantially different, as well as the efficiency of turbulent 
transport in the convective and neutron-finger unstable zones. 



Less is known about convection in rotating PNSs. The fail- 
ure of detailed simulations in spherical symmetry to get a suc- 
cessful Supernova explosion has led different authors to argue 
that multidimensional effects, such as convection or rapid rota- 
tion can be the key to the problem (see e.g. Buras et al. 2003, 
and references therein). From theoretical modeling and simple 
analytic considerations it is commonly accepted that core col- 
lapse of a rotating progenitor leads to differential rotation of 
a newly born neutron star (Zwerger & Miiller 1997, Rampp 
et al. 1998, Liu 2002, Dimmelmeier et al. 2002, Miiller et al. 
2003), mainly due to conservation of the angular momentum 
during collapse. A recent study on evolutionary sequences of 
rotating PNSs (Villain et al. 2003 1 shows that the typical scale 
on which the angular velocity changes is in the range ps 5 — 10 
km. Granted that a PNSs is born rapidly rotating, when the an- 
gular velocity is of the order or larger than the Brunt-Vaisala 
or Ledoux frequencies, even relatively fast convective insta- 
bilities can be substantially modified, and convection can be 



constrained to the polar regions (Fryer & Heger 2000 1. In this 
paper, we study this effect in the context of PNSs. 

2. The dispersion relation. 

Consider a PNS rotating with angular velocity 51 = il(s,z). 
In what follows we use cylindrical coordinates (s, ip, z). The 
deleptonization and cooling timescales are assumed to be much 
longer than the growth time of instability, thus it can be treated 
in a quasi-stationary approximation. For the sake of simplic- 
ity, we consider axisymmetric short-wavelength perturbations 
with spatial and temporal dependence exp(7i — ik ■ r) where 
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k = (k s , 0, k z ) is the wave-vector. Small perturbations will be 
indicated by a subscript 1, whilst unperturbed quantities will 
have no subscript, except when indicating vector components. 
In the unperturbed state, the PNS is in hydrostatic equilibrium, 



^=G, G = g + n 2 s 
P 



(1) 



where g is the gravity. 

In the Boussinesq approximation, the linearized momen- 
tum and continuity equations read 



71*1 + 2fl x v% + e v s(vi ■ VO) 



ikpi 



G 



Pi 



(2) 



waves in a differentially rotating fluid (see, e.g., Goldreich & 
Schubert 1967) by the presence of a term proportional to VY". 
Non-trivial solutions of Eq. @ are given by 



7 



-q 2 + u> 2 



\. (8) 

The condition of instability is u) 2 + uj 2 l > q 2 . In non-rotating 
PNSs, the characteristic frequencies in the convective regions 
are typically of the order of 1 ms _1 (MPU). Therefore, rota- 
tional effects will be important when the angular velocity £1 
is of the same order, £1 as 1000 rad/s, which can be reached 
during the early stages of PNS evolution (Villain et al.. l2003i . 



k ■ v\ = , 



(3) 3. The instability criteria. 



where V\, p\ and p\ are perturbations of the velocity, pressure 
and density, respectively; e v is the unit vector in the azimuthal 
direction. We assume that the matter is in chemical equilib- 
rium, thus the density is generally a function of the pressure 
p, temperature T and lepton fraction Y . In the Boussinesq ap- 
proximation, perturbations of the pressure are small, therefore 
pi can be expressed in terms of the perturbations of tempera- 
ture, T\, and lepton fraction, Y\, 



Pi 



i'\PY+SYi 



(4) 



where (3 and 8 are the coefficients of thermal and chemical ex- 
pansion; [3 = — {d\np/d\nT) p Y, 5 — — (<91n p/dY) p T- 

We consider instabilities arising on a timescale much 
shorter than the dissipative timescale, i.e., those perturbations 
with relatively long wavelengths but still smaller than the pres- 
sure scale, to ensure that our approximation is valid. For exam- 
ple, perturbations with A = 2n/k ~ 0.3 — 1 km satisfy both 
these requirements. For such perturbations, we can neglect dis- 
sipative effects in the equations governing thermal balance and 
lepton fraction (see MPU for more details about the influence 
of dissipative terms). Then, the linearized transport equations 
read 



Ti - Vl ■ AVT = 0, 



(5) 



Yi + vi-VY = 0, (6) 
where AVT is the super-adiabatic temperature gradient: 

( dT \ 

AVT = ( — I Vp-VT . 
\ d PJs,Y 

The dispersion equation corresponding to Eqs. is 

7 2 (7 2 + g 2 -^ 2 -ci) = o, (7) 

where 



h 2 h h 

fc 2 k 2 
= —(J3/T)G± ■ AVT 



lu 2 l = 6G±_ ■ W, 



and Gi = G — k(k ■ G)/k 2 ; fl e is the epicyclic frequency, 
Q 2 e = ^d{s 4 n 2 )/ds, and we have defined Vt 2 z = s dVL 2 /dz. 
Eq. @ differs from the standard dispersion equation of buoyant 



Assume that a vector k forms an angle <f> with the polar axis. 
Defining the following vector, 



-{/3/T)AVT + 6VY 



(9) 



the dispersion relation can be written in terms of the compo- 
nents of C and G as follows: 
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2 — A cos 2 . 



B sin d> cos < 



C Z G Z 



(10) 



where A = -n 2 e +C s G s -C z G z andB = Vl 2 -C z G s -C s G z . 
By taking the curl of Eq. 0, it can be readily obtained that the 
condition of hydrostatic equilibrium leads to 



n 2 z = [C x G] = C Z G S - C S G Z 



(11) 



so that B can be further simplified to obtain B = —2C S G Z . 

Since the dependence of Eq. (II Oi on (f> is rather simple, we 
can obtain that the maximum of 7 2 corresponds to 



cos 



1± 



A 2 



A 2 +B 2 



(12) 



with the ± sign depending on the sign of A, and the maximum 
value of 7 2 is given by 



7 2 

! max 



-n 2 e -u 2 



+ 



+- V^ 2 + toff + 4G z (c z n 2 e - c s w z ). 

where ojq — — (C ■ G). The two conditions for stability follow 
straightforwardly from the above expression: 



- o 2 - lj 2 < o ; G z (c z n z e - c s n 2 z ) < o , 



(13) 



These two conditions look like the Solberg-H0iland conditions 
(Tassoul 2000 1, but with additional terms due to the lepton 
gradients. Criteria (II 3i can be also applied to quasi-Keplerian 
toroidal configurations made of nuclear matter surrounding 
black holes, in which the effect of lepton gradients is relevant. 

Notice that if isobaric and isopicnic surfaces coincide, C 
and G are aligned, and vertical gradients of £1 are not allowed. 
Consistently, most of the studies of rotating stars with simpli- 
fied polytropic equations of state use angular rotation profiles 
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Fig. 1. Characteristic squared frequency normalized to the ro- 
tation frequency, 7 2 /f2 2 , as a function of the angle between 
the wavevector and the polar axis (0), assuming that uJq = fi 2 . 
The rotation law is described by Eq. dl4> . Three different radial 
points, with s/Rq = 0.1, l,and 5 are plotted with solid lines, 
dots, and dashes, respectively. 

that depend only on s. In our calculations, we use the follow- 



ing law proposed by Komatsu et al. {1989 1 and used by many 
authors: 



fi(s) 



(14) 



In Fig. 1 we show j 2 /fl 2 as a function of the wavevec- 
tor orientation in a convectively unstable situation, ojq > 0. We 
have taken ujq — £! 2 , and the polar angle 9 = n/4. For simplic- 
ity, we have used the law (114-1 for the angular velocity so that 
Vl z = 0. Three different cases are plotted with solid lines, dots, 
and dashes, corresponding to s/Rq = 0.1,2, and 5, respec- 
tively. When s /Rq <C 1 the angular velocity varies slowly and 
the situation is close to rigid rotation, while for s /Ro ^> 1 the 
specific angular momentum is nearly constant and rotational 
effects are negligible. When the angular momentum increases 
with the distance to the axis (Rayleigh stable), j 2 can become 
negative and instability is suppressed for a subset of the possi- 
ble orientations of k. Only perturbations with k oriented more 
or less perpendicular to the polar axis (fluid motions in the ver- 
tical direction) are not suppressed. Note that the fact that ro- 
tation results in turbulent convective motions preferentially in 
the direction parallel to the rotation axis was first pointed out by 
Randers ( 1942 > in the context of the study of energy transport 
in stellar convective cores. 

This effect is more manifest looking at Figs. |2] and [5] In 
Fig. 2 we plot the maximum value of 7 2 as a function of 
the polar angle (9), at three different spherical radii. Near the 
pole the maximum of 7 does not change, but as we approach 
the equator, rotational effects may reduce its value drastically, 
unless the gradient of the angular momentum is small (case 



Fig. 2. Maximum of 7 2 /SI 2 as a function of the polar angle for a 
fixed radius of r/Ro = 0.1 (solid), 2.0 (dots) and 5.0 (dashes). 

t/Rq = 5). More importantly, not only does the growth time 
become longer, but the unstable perturbations are restricted to 
some range of angles, centered at the value given by Eq. H2\ . 
and width given by 



tan Ac/) 



■ w 



(15) 



In Fig.|5]we plot surfaces of constant 7 2 /Sl 2 in a 9—<f) plane, for 
r/Ro = 1. The curve 7 2 = encloses the interval of angles in 
which perturbations are unstable. The position of the maximum 
is indicated by the dashed line. Notice that when the second 
of the stability conditions (II 31 is not fulfilled, the function 7 2 
has no zeros. Therefore, 7 2 is either positive or negative for all 
angles depending on the sign of the first condition. 

4. Discussion. 

We have derived the criteria for convective instability in rotat- 
ing, non-magnetic PNSs, obtaining two conditions analogous 
to the Solberg-H0iland criteria but including the effects of lep- 
ton gradients. In convectively unstable regions, the existence of 
a positive gradient of angular momentum can stabilize the un- 
stable convective modes in some directions. The basic picture 
is that an unstable buoyant element of fluid is stabilized when it 
tries to move toward a direction in which the surrounding fluid 
elements have larger angular momentum, because it is subject 
to a smaller centrifugal force and therefore restored to its equi- 
librium position. However, this restoring force has no effect on 
buoyant waves that follow lines of constant angular momen- 
tum. For example, in a situation in which thermal, lepton, and 
pressure gradients are nearly aligned, the angular momentum 
depends weakly on the vertical coordinate, thus allowing con- 
vectively unstable modes to grow only in the direction approx- 
imately parallel to the rotation axis. In a general case, the ge- 
ometry will be more complicated but the qualitative results will 
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Fig. 3. Surfaces of constant j 2 / ^l 2 in a 8 — <p plane, for t/Rq = 
1. The region enclosed by the surface 7 2 /0 2 = is convec- 
tively unstable. The dashed line indicates the angle at which 
the maximum is reached. 

not change. In Fig.|4]we illustrate this effect in a realistic model 
of a PNS 0.5 s after birth, taken from Pons et al. (1999). The 
convectively unstable region extends from 17 to 28 km, and the 
arrows indicate the direction of the most rapidly growing mode, 
with the length being proportional to the inverse growth time. 
In the convectively unstable region near the equator, fl e varies 
between 500 and 1000 rad/s, which is about one third the value 
of uiq, and we have taken Ro = 10 km (defined in Eq. H41 ). In 
a saturated state, hydrodynamical non-linear effects will proba- 
bly smear out anisotropies on small scales but it is unlikely that 
this may occur on the largest convective scale. Hydrodynamic 
simulations of the early evolution of rotating PNSs show that 
the scenario drawn in Fig.|4]is fulfilled during the very early 
phase of convection in PNSs (H-Th. Janka, private communi- 
cation), which shows that our linear stability analysis reflects 
accurately the situation at the onset of convection. However, af- 
ter several hundreds of milliseconds, the efficient neutrino and 
convective transport in the polar region will flatten the temper- 
ature, lepton and angular momentum profiles, resulting in a dif- 
ferent situation: convection occurring only close to the equator, 
although convective shells are still elongated and parallel to the 
rotation axis (see e.g. Fig. 10 in Janka & Keil ll997l or Janka, 
Kifonidis & Rampp [T997h . 

This coupling between convective instability and stabi- 
lizing gradients of angular momentum can be another in- 
gredient that helps to explain the large asymmetry that has 
been observed in recent simulations of rotational core collapse 
Supernovae (Miiller et al., 2003), not observed in simulations 
of convective PNSs without rotation. The mechanism we pro- 
pose is a natural way to create anisotropic energy and momen- 
tum transport by convective motions, that only requires that the 
angular velocity be of the order of the Brunt-Vaisala or lep- 
tonic frequencies. Notice that our results are relevant for con- 
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Fig. 4. The direction of the motion for the most rapidly grow- 
ing modes in a realistic model of PNS 0.5 s after birth. In the 
convective region (17-28 km), the instability is more effective 
in the polar region and, even at intermediate latitudes, vertical 
motions are preferred. 

vective motions in the PNS, more than in the extended enve- 
lope as in Fryer & Heger (2000 1. Our analysis is only valid 
in a steady-state situation, while the conditions in the accre- 
tion layer are extremely non-stationary, and other effects such 
as anisotropic accretion are certainly more relevant. If the col- 
lapsed core shows strong differential rotation, at scales larger 
than 50-100 km the angular momentum can be nearly constant 
and rotational effects are small. Nevertheless, the anisotropy 
originated in the inner core may have some effect at larger 
scales, and can help to explain the Supernova mechanism. 
Obviously, self-consistent 3D simulations involving convec- 
tion, neutrino transport and angular momentum transport are 
needed prior to drawing robust conclusions. 

In our study we only considered non-magnetic proto- 
neutron stars. The presence of a large magnetic field could 
change this picture. It has recently been shown that turbulent 
mean-field dynamo can be effective for PNSs with periods 
shorter than 1 s, which would generate strong magnetic fields in 
the interior (Bonanno et al. 2003). However, magnetic field ef- 
fects begin to be relevant only when the Alfven frequency is of 
the order of the other dynamical or dissipative frequencies, and 
this requires magnetic fields larger than 10 15 G (Miralles et al. 
2002 1. Probably the missing ingredient to explain Supernovae 
is related to the coupling between different effects, such as the 
one proposed in this paper, the magneto-rotational instability 
(Akiyama et al., 2003}, or a combination of all of them. 
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